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New  Procedure  for  Computing  Finite-Amplitude  Distortion 

Bill  D.  Cook 

Dtpottm»nt  of  Physics,  MicUgon  Slate  University,  East  Lansing,  Michigan 
(Received  November  21,  1961) 

An  iterative  proceas  is  described  for  the  calculation  of  the  distortion  of  plane  6nite-amplitude  sound 
waves  in  a  dissipative,  nondisperuve  medium.  The  method  of  calculation  is  a  discrete-interval  process  of 
considering  the  distortion  of  the  wave  while  propagating  through  a  small  distance,  correcting  for  absorption 
within  this  distance,  and  then  considering  tUs  new  wave,  etc.  It  is  necessary  to  use  a  high-speed  electronic 
computer  to  obtain  the  spectral  composition  of  the  wave.  The  iterative  process  allows  calculations  beyond 
the  “discontinuity  distance.”  The  spatial  change  of  the  spectra]  composition  is  used  for  the  calculation  of 
“absorption  coefficients”  describing  the  energy  loss  from  the  total  wave  and  from  the  fundamental  com¬ 
ponent.  These  absorption  coefficients,  which  are  functions  of  distance  from  the  source,  arc  found  to  be 
remarkably  different. 


INTRODUCTION 

N  recent  years  the  problem  of  finite-amplitude  sound 
waves  has  been  approached  through  the  solution  of 
the  nonlinear  equation  of  motion.  Although  this  ap¬ 
proach  has  answered  some  of  the  questions  concerning 
the  propagation  of  finite-amplitude  waves,  the  range 
of  reliability  of  these  theories,  which  are  usually  based 
on  equations  approximating  the  exact  differential  equa¬ 
tion,  is  small.  The  reliability  of  this  usually  extends  up 
to  the  neighborhood  the  "discontinuity  distance.” 
However,  from  the  p.ienomenological  aspects  of  the 
problem,  one  can  formulate  a  very  simple  model  from 
which  solutions  at  greater  distances  may  be  obtained. 
Presented  here  is  such  a  model.  While  it  necessitates 
the  use  of  a  high-speed  electronic  computer,  it  allows 
one  to  compute,  by  an  iterative  process,  the  spectral 
composition  of  the  waves  at  all  distances.  It  is  not  the 
purpose  of  this  paper  to  correlate  the  results  of  previous 
theoretical  and  experimental  investigations,  but  rather 
to  show  what  may  be  obtained  from  a  very  simple 
model.  The  readers  are  referred  to  the  papers  of  Keck 
and  Beyer*  and  Zarembo  and  Krasil’nikov’  for  reviews 
of  existing  analyses  using  the  differential  equation. 

A  model  approach  was  previously  used  by  Fox  and 
Wallace,'  who  used  a  graphical  analysis  to  determine 
the  spectral  composition  and  consequently  determined 
an  “absorption  coefficient”  for  finite-amplitude  waves. 

■  W.  Keck  and  R.  T.  Beyer,  Phys.  Fluids  3,  346-352  (1960). 

’  L.  K.  Zarembo  and  V.  A.  Krasil'nikov,  Soviet  Phys. — Uspekhi 
2,  580-599  (1959). 

'  F.  E.  Fox  and  W.  A.  Wallace,  J.  Acoust.  Soc.  Am.  26, 994-1006 
(1954). 


Although  the  model  here  is  based  somewhat  on  the 
same  assumptions  as  those  of  the  Fox  and  Wallace 
theory,  a  more  basic  aspect  of  the  absorption  mecha¬ 
nism  is  introduced.  A  high-speed  electronic  computer  is 
used  to  calculate  the  distortion.  A  parametric  integra¬ 
tion  method  similar  to  that  given  by  Fubini-Ghiron' 
is  used.  This  integral  method  of  Fubini-Ghiron  has 
also  been  given  by  Hargrove'  and  Keck  and  Beyer.' 

This  model  is  based  on  two  fundamental  assumptions 
which  describe  the  distortion  and  absorption  mecha¬ 
nisms.  These  assumptions  are  the  following ; 

(1)  The  distortion  mechanism  can  be  described  by 
a  change  in  phase  velocity  directly  proportional  to  the 
particle  velocity. 

(2)  The  absorption  mechanism  can  be  described  by 
assuming  that  the  rate  of  the  absorption  of  each 
harmonic  is  proportional  to  the  amount  of  the  harmonic 
present  and  to  the  square  of  the  frequency  of  the 
harmonic. 

These  basic  postulates  are  applied  to  a  wave  by  first 
allowing  the  wave  to  distort  while  it  propagates  through 
a  small  interval,  and  then  correcting  for  absorption. 
This  wave  of  new  shape  is  now  allowed  to  distort  and 
be  absorbed.  By  assuming  that  no  discontinuity  of  the 
wa're  shape  is  formed  because  of  absorption,  one  may 
calculate  by  this  continuing  process  the  shape  of  the 
wave  at  ail  distances.  Although  the  absorption  and 
generation  of  the  harmonics  are  treated  indepiendently 


*  R.  Fubini-Ghiron,  Alta  Frequenz.  4,  530  (1935). 

>  L.  E.  Hargrave,  J.  Acoust.  Soc.  Am.  32,  511-512  (1960). 
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Fig.  1.  Diagram  showing  the  distortion  process 
postulated  for  a  finite-amplitude  wave. 


in  each  small  interval,  the  end  result  contains  the  inter¬ 
action  between  absorption  and  generation  mechanisms. 

Others*'^  have  considered  similar  absorption  mecha¬ 
nisms  along  with  distortion ;  however,  the  effect  of  the 
absorption  on  the  generation  of  the  harmonics  is  neg¬ 
lected.  The  range  of  these  theories  which  are  based  on 
the  Fubini-Ghiron  method  does  not  exceed  the  dis¬ 
continuity  distance  as  given  by  the  dissipationless 
theory. 

The  model  presented  here  is  not  expected  to  be  valid 
for  extremely  intense  sound  waves  which  may  produce 
periodic  shock  fronts  but  rather  for  sound  waves  of 
moderate  amplitudes  traveling  in  a  fairly  absorbing 
medium.  In  particular,  this  model  .should  be  useful  for 
the  investigation  of  finite-amplitude  waves  of  the  mega¬ 
cycle  frequencies  in  liquids.  The  calculations  presented 
from  this  model  are  for  a  range  of  parameters  which 
cover  the  practical  experimental  situations. 


FORMULATION  OF  THE  MODEL 

In  this  paper  we  shall  consider  plane  acoustic  waves 
in  a  nonlinear,  nondispersive  medium.  Following  the 
usual  notation  we  write  an  equation  relating  the  in¬ 
stantaneous  pressure  p  and  density  p  as 


P=P<\-\-A 


The  zero  subscripts  refer  to  the  undisturbed  medium. 
The  terms  A  and  B  are  constants  for  a  given  medium 
at  a  given  temperature.  Using  the  above  equation  of 
state,  the  phase  velocity  is 

c=co-f[B/(2^)-|-l]«,  (2) 

accurate  to  the  first  order  in  u,  where  u  is  the  particle 
velocity  and  co  is  the  velocity  of  sound  with  infinitesimal 
amplitude.  It  is  assumed  here  that  this  linear  change 
of  phase  velocity  with  particle  velocity  is  the  only 
mechanism  which  causes  the  wave  to  distort.  Points  on 
the  wave  having  larger  values  of  particle  velocity  tend 
to  overtake  points  of  lesser  values.  For  a  dissipationless 
medium,  an  initially  sinusoidal  waveform  will  become 
multiple  valued  at  distances  greater  than  a  distance  L 

•  W.  W.  1-ester,  J.  Acoust.  Soc.  Am.  33,  1196  (1961). 

’L.  Adler,  "A  Study  of  the  Distortion  of  Finite  Amplitude 
Ultrasonic  Waves  in  Liquids,”  M.S.  thesis,  Michigan  State  Uni¬ 
versity,  East  Lansing,  Michigan,  1961. 


given  by* 

L=c^l2[B/(2A)+l]/r[/:(0))-',  (3) 

where  /  is  the  frequency  of  the  wave  and  Ut(0)  is  the 
maximum  particle  velocity  of  the  initial  wave.  This 
distance  L  (known  as  the  “discontinuity  distance”) 
shall  become  one  of  the  important  parameters  in  de¬ 
scribing  the  results  from  this  model  approach. 

The  exact  frequency  dependence  of  the  absorption 
coefficient  for  most  media  is  unknown.*  Here  we  shall 
assume  for  simplicity  that  the  rate  of  absorption  is 
proportional  to  the  square  of  the  frequency  of  the 
wave ;  this  is  characteristic  of  a  viscous,  heat-conducting 
medium. 

In  the  absorption  of  infinitesimal  waves,  the  expo¬ 
nential  decay  which  describes  the  “diffusive  nature” 
of  the  absorption,  occurs  because  the  damping  forces 
produced  are  proportional  to  the  quantity  being 
damped.  Although  for  acoustic  waves  in  a  viscous, 
heat-conducting  medium,  the  exponential  decay  arises 
from  a  higher-order  differential  equation,  we  shall  as¬ 
sume  that  the  “diffusive  nature”  can  be  described  by 
a  first-order  differential  equation,  namely,  the  rate  of  ab¬ 
sorption  of  a  quantity  is  proportional  to  that  quantity. 
However,  in  considering  the  “diffusive”  absorption  of 
waves,  the  exponential  is  a  special  case.  It  requires 
that  the  quantity  being  absorbed  will  remain  at  a  con¬ 
stant  value  if  the  absorption  mechanism  vanishes. 
However,  since  the  amplitudes  of  the  harmonics  are 
changing  due  to  distortion  while  undergoing  absorption, 
the  diffusive  absorption  process  must  be  considered  for 
each  harmonic  according  to  its  growth.  The  harmonics 
are  treated  independently,  i.e.,  the  absorption  of  a 
harmonic  in  a  wave  is  the  same  as  the  absorption  of  an 
infinitisimal  wave  of  that  harmonic.  Details  of  how  the 
generation  affects  the  rate  of  the  absorption  are  given 
in  the  Appendix. 

With  these  distortion  and  absorption  mechanisms, 
it  is  relatively  easy  to  calculate  the  spectral  composi¬ 
tion  of  a  finite  amplitude  wave  by  an  iterative  process. 
In  such  a  process,  the  results  are  expected  to  be  more 
valid  for  smaller  increments  of  interval.  No  anal)rtical 
treatment  will  be  given  here  for  the  estimate  of  size 
of  these  intervals  in  terms  of  the  other  parameters.  It 
should  be  noted,  however,  that  for  small  absorption 
parameters  or  greater  amplitudes  of  the  wave,  the 
wave  approaches  more  closely  to  the  dissipationless 
case.  For  these  conditions,  it  is  reasonable  to  assume 
that  intervals  must  be  small. 

DETAILS  OF  COMPUTATIONS 

To  calculate  the  effects  of  distortion,  we  shall  now 
gfneralize  the  method  outlined  by  Hargrove*  for  the 
dissipationless  case.  Let  a  function  u=f{x)  describe 
the  particle  velocity  at  a  given  time.  The  independent 

'  J.  J.  Markham,  R.  T.  Beyer,  and  R.  B.  Lindsay,  Revs.  Modem 
Phys.  23,  353-411  (1952). 
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variable  is  a  reduced  spatial  coordinate  which  describes 
the  function  in  the  range  0<*<2ir.  For  simplicity  we 
shall  now  restrict  /(*)  to  be  an  odd  function.  For 
arbitrary  waveforms  this  restriction  can  be  removed. 
During  a  small  interval  of  time  At,  the  whole  waveform 
propagates  a  distance  CoAI.  However,  in  addition  to 
this  velocity  co,  each  point  of  the  wave  described  by 
the  function  u  is  assumed  to  propagate  at  a  velocity 
proportional  to  the  value  of  the  function.  The  velocity 
of  each  point  is  given  by  Eq.  (2).  Figure  1  shows  one 
half  of  an  arbitrary  waveform  [described  by  /(*)]  and 
the  resultant  waveform  [described  by  /i(*)]  from 
such  a  distortion  mechanism.  The  problem  now  is  to 
calculate  the  function  fi{x)  provided  that  f(x)  and 
B/A  are  known. 

As  f(x)  is  an  odd  function,  fi(x)  is  necessarily  an 
odd  function.  Thus  both  can  be  expanded  in  a  Fourier 
sine  series.  Let 

fix)^^^  b„smnx  (4) 

ir-l 

and 

oe 

/iW=L  6»,tsin«x.  (5) 

If  fi{x)  is  known,  the  expansion  coefficients  can  be 
calculated  from 

2  r 

/  fi(x)smttxdx.  (6) 

r  Jo 

Graphically,  one  may  obtain  fi{x)  if  f(x)  is  known  and 
thus  obtain  the  coefficients  6„,i.  By  parametric  inte¬ 
gration,  one  may,  in  principle,  analytically  find  these 
coefficients  if  /(*)  and  B/A  are  known. 

To  evaluate  Eq.  (6)  let 

«=/[*(»)],  (7) 

WnPrP 

x{e)=e+km-  (8) 

The  term  kj{8)  describes  the  distortion  mechanism 

as  required  by  Eq.  (2).  The  value  of  k  is  restricted  to 
be  sufficiently  small,  such  that  the  waveform  does  not 
become  multiple  valued.  Hargrove  has  shown  that  k 
represents  the  fractional  distance  to  the  discontinuity 
distance  L  of  an  initially  sinusoidal  wave.  In  other 
words,  if  in  a  dissipationless  medium  an  initially 
sinusoidal  wave  propagates  through  a  distance  R,  the 
spectral  composition  can  be  calculated  using  k=R/L. 
Substitution  of  Eqs.  (7)  and  (8)  into  Eq.  (6)  gives 

bn.i=- f  {miinlne+nkf(emi+kf'{e))de.  (9) 

It  Jo 

Thus  if  we  know  the  description  of  the  particle  velocity 
at  any  given  time,  we  can  find  the  description  of  the 
particle  velocity  at  a  new  time  if  we  are  considering 
only  the  distortion  mechanism. 


Fig.  2.  Normalized  spec¬ 
tral  composition  of  a  typical 
finite-amplitude  wave  versus 
the  reduced  distance  K. 


The  procedure  for  introducing  absorption  is  the 
following ;  assume  a  waveform  (described  by  the  coeffi¬ 
cients  choose  a  small  increment  of  k\  calculate 
the  spectral  com]x>sition  of  a  distorted  wave  by  Eq. 
(9)  (in  terms  of  coefficients  correct  the  coeffi¬ 

cients  (then  described  by  6n,<+i')  for  absorption;  and 
repeat  the  above  process  using  the  previous  results  to 
describe  the  new  wave.  Each  coefficient  is  cor¬ 
rected  independently  according  to  the  assumption  that 
the  absorption  throughout  the  small  interval  of  k  is 
proportional  to  the  amount  of  the  harmonic  present 
at  any  given  point  in  the  interval.  The  frequency 
dependence  of  the  absorption  must  be  considered. 

From  the  dissipationless  theory, 

ft,=  (-l)»(2/»A)/,(«/fe).  (10) 

It  can  be  seen  that  A6,,((=6„,<+i— is  approximately 
described  by  the  first  term  of  the  power-series  expres¬ 
sion  of  Eq.  (10)  in  the  f  orm  Aft  for  small*. 

Assuming  this  power  dependence  on  k,  the  absorption- 
corrected  coefficients  can  be  written  as 

b„.i^i'-bn.i!tx^{-aon'^Lk)+Abn.Af>-i{aon^Lk).  (11) 

The  first  term  of  the  right-hand  side  of  Eq.  (11)  is  the 
exponential  decay  if  Aft„,i=0;  oo  is  the  absorption 
coefficient  of  a  wave  of  the  fundamental  frequency  and 
infinitesimal  amplitude.  The  function  .4»_i(ao«®L*)  cor¬ 
rects  for  the  absorption  of  the  changing  function  Aftn.i. 
This  function  A  »_i,  which  is  described  in  the  Appendix, 
is  a  pwwer  series  for  which  the  exponential  power  series 
is  a  special  case.  The  »-squared  dependence  of  the  ex¬ 
ponential  and  An-i  is  introduced  by  the  assumption 
of  the  frequency  squared  dependence  of  the  absorp¬ 
tion.  The  product  kL  gives  the  actual  distance  that  the 
wave  propagates  for  the  increment  k. 

For  the  calculations  presented  here  the  initial  wave¬ 
form  is  assumed  to  be  sinusoidal  with  the  unit  ampli¬ 
tude  (fti,o=  1).  The  corresponding  particle  velocity  and 
pressure  may  be  determined  through  the  parameter  L. 
The  results  are  expressed  in  terms  of  the  product  aoL 
and  the  reduced  distance  K.  This  distance  K,  which 
is  the  distance  that  the  wave  has  prop>agated  from  the 
origin,  is  given  by  K=mk,  where  m  is  the  number  of 
intervals.  When  ^r=l,  the  wave  has  propagated  a 
distance  equal  to  the  discontinuity  distance. 
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Fia.  3.  Companion  of  the 
normalized  inteiMity(ioUd  line) 
calculated  from  spectral  curves 
of  Fig.  2  and  the  exponential 
decay  expected  for  a  wave 
undergoing  no  finite-amplitude 
distortion. 


where  I  is  the  intensity  of  the  sound  wave.  This  coeffi¬ 
cient  (XT  describes  the  rate  at  which  the  energy  of  the 
total  wave  is  dissipated  at  any  given  distance.  One 
may  also  define  an  absorption  coefficient  which  de¬ 
scribes  the  rate  of  loss  of  energy  from  the  fundamental 
comp>onent  alone  as 

a,=  -rf(f7,*)/2t' -dUi/U4x. 

Figure  4  shows  these  coefficients  (normalized  to  ao)  for 
aoL=0.10.  It  is  important  to  stress  that  the  absorption 
coefficient  for  the  fundamental  frequency  ai  differs 
from  the  total  absorption  coefficient  ay  for  most  values 
of  K.  It  can,  therefore,  be  concluded  that  the  total 


The  calculations  were  done  by  the  Michigan  State 
integral  computer.  To  accurately  retain  the  shape  of 
the  waveform,  as  many  as  16  harmonics  were  computed 
for  each  interval.  In  the  results  presented  here,  the 
value  of  k  never  exceeds  O.OS,  one-twentieth  of  the  dis¬ 
continuity  distance. 

DISCUSSION  OF  RESULTS 

Figure  2  shows  the  spectral  composition  of  a  wave 
from  a  typical  calculation  using  this  model.  This  com¬ 
position  is  given  in  terms  of  the  harmonic  structure  of 
the  particle  velocity.  The  amplitudes  Un{K)  of  the 
harmonics  are  normalized  to  the  particle  velocity  am¬ 
plitude  f7i(0)  of  the  initial  sinusoidal  wave.  The  value 
of  aoZ,  equals  0.10  for  these  calculations  which  extend 
to  10  times  the  discontinuity  distance.  For  this  cal- 


Fig.  4.  Normalized 
absorption  coeffici¬ 
ents  for  the  total 
wave  nr/a«  and  for 
the  fundamental 
component  ai/ag  for 
the  values  given  in 
Figs.  2  and  3. 


culation,  the  amplitudes  of  higher  harmonics  peak  in 
the  region  1  <K  <2. 

Figure  3  shows  the  effect  of  the  nonlinear  propaga¬ 
tion  on  the  dissipation  of  the  sound  intensity  calculated 
from  the  curves  presented  in  Fig.  2.  As  is  expected, 
there  is  a  large  difference  between  an  exponential  decay 
and  the  decay  predicted  by  this  model.  The  predicted 
greater  total  absorption  is,  of  course,  caused  by  the 
higher  rate  of  absorption  of  the  generated  harmonics. 
To  describe  this  higher  rate  of  absorption,  one  may 
define  a  finite  amplitude  absorption  coefficient'  as 


Fig.  S.  Normalized 
total  absorption  coeffi¬ 
cients  for  various  values 
of  aoL. 


absorption  coefficient  is  not  determined  by  expierimental 
absorption  measurements  of  the  fundamental  compo¬ 
nent  only.  The  difference  between  coefficients  oti  and 
oir  can  be  very  large,  especially  at  small  distances; 
the  energy  of  the  fundamental  is  being  lost  by  both 
distortion  and  absorption  while  the  energy  of  the  total 
wave  is  being  lost  only  by  absorption. 

As  to  be  expected,  the  values  of  ar/ao  and  oi/oio 
depend  on  the  parameter  ooL.  Figures  5  and  6  show 
the  results  for  or/oo  and  oi/ao,  respectively,  for  various 
values  of  the  parameter  ooL.  For  large  values  of  aoZ. 
(larger  than  0.50),  which  means  either  high  absorption, 
small  pressure,  or  small  nonlinearity,  the  absorption 


Fio.  6.  Normalized 
fundamental  abioiption 
coeffident  for  various 
values  of  atL. 
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tend  to  become  negligible  even  at  moderate  values  of 

a^L. 

CONCLUSION 

An  iterative  process  is  used  to  determine  the  spectral 
composition  of  plane  finite-amplitude  waves  in  a 
dissipative  medium.  From  the  spatial  change  of  the 
amplitudes  of  the  harmonic  components,  “absorption 
coefiScients”  describing  the  loss  of  the  energy  from 
the  total  wave  and  from  the  fundamental  component 
are  computed.  The  results  are  presented  as  families 


Fio.  7.  The  maximum  values  of  the  oonnalixed  aheorption 
coeffidents  versus  the  parameter  cttL. 

coefficients  ar  and  at  are  nearly  equal  to  as.  However, 
for  small  values  of  atL,  these  absorption  coefficients 
become  large  compared  to  as.  Also  at  small  values  of 
K,  the  absorption  coefficients,  ar  and  at,  differ  greatly. 
For  example,  for  Ar=O.SO  and  asL=0.01,  the  ratio 
ai/ar  is  approximately  equal  to  10. 

The  maximum  values  of  ar/as  and  a^/at  for  values 
of  aoL<0.10  occur  in  the  neighborhood  of  X=2.0,  and 
iir=l.S,  respectively.  The  dependence  of  the  maximum 
values  of  ar/ao  and  ai/as  on  aoL  is  shown  in  Fig.  7. 

In  Fig.  8  the  fundamental  component  is  shown  for 
several  values  of  asZ,.  As  atL  decreases,  the  values  of 
Ui(K)/Ui(())  are  approaching  upper  limits.  Of  course. 


Fig.  8.  Normalized 
fundamental  compo¬ 
nent  for  various 
values  of  ooL.  Dashed 
line  gives  values 
of  dissipationless 
theory. 


of  curves  having  the  parameter  atL.  This  parameter 
determines  the  experimental  situation  as  ao,  the  ab¬ 
sorption  coefficient  of  the  fundamental  component 
of  infinitesimal  amplitude,  and  L,  the  “discontinuity 
distance,”  describes  the  absorption  and  distortion 
mechanisms. 
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for  K<1  the  limiting  values  are  those  given  by  the 
dissipationless  theory,  not  only  for  the  fundamental, 
but  also  for  the  higher  harmonics. 

Figure  9  shows  the  effect  of  the  parameter  atL  on 
the  second  harmonic.  For  values  of  aoL<0.20,  the 
maximum  value  of  the  second  harmonic  occurs  at 
K<»l.2S.  For  larger  values  of  atL,  these  maxima  are 
seen  to  occur  at  lower  values  of  K.  Similarly,  the 
maximum  values  of  the  third  and  fourth  harmonics 
are  affected  by  the  parameter  atL  as  shown  in  Figs. 
10  and  11.  As  a  consequence  of  the  frequency  squared 
dependence  of  the  absorption,  these  higher  harmonics 


Fio.  10.  Normal- 
ized  third  harmonic 
for  varioui  values  of 
atL.  Daihed  line 
gives  values  of  dis- 
npationless  theory. 
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Fig.  11.  Normalized 
fourth  harmonic  for 
various  values  of  atL. 
Dashed  line  gives  values 
of  dissipationless  theory. 


The  work  was  suppiorted  by  the  U.  S.  Army  Research 
Office  (Durham)  and  by  the  Office  of  Naval  Research. 

APPENDIX 

By  assuming  the  rate  of  loss  of  a  quantity  to  be 
proportional  to  the  quantity  present  at  each  instant, 
one  expects  a  decay  of  the  quantity  according  to  an 
exponential  law.  However,  this  is  true  only  if  the 
quantity  would  remain  constant  if  the  loss  mechanism 
were  absent.  We  shall  consider  here  the  case  where  the 
quantity  changes  according  to  a  power  law  in  the 
interval  if  absorption  were  not  present. 

Let  the  dependence  of  the  quantity  y  on  the  inde¬ 
pendent  variable  *  be  y=Cx’'  if  there  were  no  absorp¬ 


tion.  By  assuming  the  loss  mechanism  described  above, 
one  can  write  the  integral  equation 


(Al) 


where  a  is  a  constant.  This  integral  equation  is  of  the 
type  known  as  Volterra’s  linear  integral  equation  of 
the  second  kind;  the  solution  may  be  obtained  by 
successive  substitutions.'  The  solution  obtained  by  this 
method  may  be  written  as 

y=C*”Am(o*),  (A2) 

where  Amiax)  is  the  hypergeometric  series 


ax  (axY 

^«(o*)  =  l - 1 - 

»n-|-l  (»»-|-l)(f»-t-2) 

(w-f-1)  (♦»-(- 2)  (w-t-3) 


(A3) 


This  series  becomes  the  exponential  series  when  >«=0, 
i.e.,  when  y  is  a  constant  if  the  absorption  were  absent. 


'W.  Lovitt,  Linear  Inleeral  Eiuations  (Dover  Publications, 
Inc.,  New  York,  19S0),  p.  13. 
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Abatract:  Iha  diatortion  of  plana  finita  anplitude  aound  wavaa  in  a  diaaipativa,  non- 
diaparaiva  aaditoi  can  ha  calculated  by  an  itarativa  procaaa  of  allowing  tha  wave  to 
dlatort  and  than  corracting  for  tha  abaorptlon.  gaaulta  of  auch  calculationa  ara 
praaantad. 

In  tha  laat  few  yaara,  thara  haa  baan  conaidarabla  Interaat  in  tha  problaa  of 
flnlta-aaplitude  aound  wavaa  in  diaaipativa  aMdla,  primarily,  bacauaa  of  tha  exceaa 
abaorptlon  of  thaae  wavaa.  Moat  theoretical  attaapta  to  calculate  finite-amplitude 
diatortion  and,  conaaquantly,  the  abaorptlon  have  not  been  entirely  auccaaaful.  The 
chief  difficulty  with  thaaa  theoretical  attampta,  which  are  uaually  baaed  on  equationa 
approximating  tha  exact  differential  equationa,  la  that  the  range  of  reliability 
uaually  extenda  only  up  to  the  neighborhood  of  the  dlacontlnuity  dlatance  aa  given 
by  the  diaalpatlonleaa  theory.  However,  from  tha  phenomenological  aapacta  of  the 
problem,  a  aimple  modal  can  be  formulated  from  which  aolutiona  for  greater  diatancaa 
can  be  calculated.  In  thia  paper,  an  outline  of  thla  model  ia  given  and  aoma  of  the 
more  Intereatlng-  concluaiona  are  noted. 

Thla  model  baaically  conalata  of  applying  diatortion  and  abaorptlon  aNchanlama 
in  an  iterative  manner.  Thia  modal  la  baaed  on  two  fundamental  proceaaea  which 
deacrlbe  the  diatortion  and  abaorptlon  mechaniama  independently.  The  mechanlam 
for  the  diatortion  la  deacribed  by  a  change  in  phaae  velocity  directly  proportional 
to  the  particle  velocity.  Thia  way  be  written  aa 


c^  +  [(B/EA)  +  1]  u 


(1) 


where  c  la  the  velocity  of  aound  wavea  of  Inflniteaimal  amplitude  and  u  la  the 
InatantSnaoua  particle  velocity.  The  tenaa  A  and  B,  which  are  conatanta  for  a  given 
medium  and  a  given  teaperature,  are  defined  by  the  equation  of  atate  relating  the 
Inatantaneoua  preaaure  p  and  denaity  p  in  the  following  equation: 


+  A  fP'Po  +  (B/2)  P  'Pq 

I  Po  J  I  Po  . 


(2) 


where  the  aero  aubacripta  refer  to  the  undiaturbed  medium.  For  thla  model  it  la 
aaauamd  that  only  thia  linear  change  in  phaae  velocity  cauaea  the  wave  to  dlatort. 
Thia  diatortion  la  produced  aa  polnta  on  the  wave  having  larger  valuea  of  particle 
velocity  tend  to  overtake  pointa  of  leaaar  valuea.  (aee  Fig.  1). 

The  abaorptlon  machaniam  la  deacribed  by  aaauming  that  the  rate  of  abaorptlon 
may  be  obtained  by  calculating  the  abaorptlon  of  each  harmonic  independently  accord¬ 
ing  to  the  amount  of  tha  harmonic  preaant  and  to 
the  aquara  of  the  frequency  of  the  harmonic.  Thia 
aaawaaa  that  the  abaorptlon  mechaniaaM  are  linear. 
With  thaaa  two  mechaniama,  tha  procedure  for  the 
calculation  of  finita-amplltude  diatortion  ia  the 
following:  aaauam  an  initial  waveform  f(x)  which 
may  be  deacribed  by  a  aat  of  Fourier  coefficienta 
b  ;  allow  tha  wave  to  dlatort  while  propagating  a 
email  dlatance,  producing  the  function  fi(x) 

deacribed  by  the  coefficienta  b  correcting 

n,  1 
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coaffleionta  for  abaorptlon;  and  than  allowing  tbaaa  eoafflelanta  to  daacrlba  tba 
now  wavafora  and  eontintlng  thia  proeaaa  In  an  Itaratlwa  wumar.  By  aaaualng  that 
tba  wavaahapa  doaa  not  bacoM  aultlpla  valuad  bacauaa  of  tba  abaorptlon,  ona  way  eal- 
culata  by  thia  continuing  proeaaa  tba  ahapa  of  tba  wava  at  all  diatancaa.  Although 
tba  ganaratlon  and  abaorptlon  of  harwonlca  aro  traatod  Indapandantly  In  aach  Intarval, 
tba  and  raault  eontalna  tba  Intaractlon  batwaan  tbaaa  wachanlaaM. 

A  wodal  aoMwhat  baaad  on  tba  aaaa  aapaeta  outllnad  hara  waa  uaad  by  Fox  and 

Wallaca(l)  in  thair  attaa^t  of  calculating  abaorptlon  eoafflelanta  for  flnito>awpll- 

tuda  wavaa.  Thay  graphically  analyaad  tba  wavaforaa  to  obtain  tba  Fourlar  eoafflelanta 

Howavar,  It  la  poaalbla  to  uaa  a  paraaatrle  Intagratlon  aathod  for  calculating  tba 

dlatortlon  proeaaa  which  allowa  tba  problaa  to  ba  anally  adaptad  to  high  apaad  cow- 

putatlon.  If  tba  function  f(x)  la  known  analytically,  tba  function  f|(6)  way  bo 

obtainad  by  Introducing  x  >  6  -f  k  f(d)  Into  f(x).  Iha  tana  kf(0)  Intfodueaa  tba 

dlatortlon  wachanian  ragulrad  by  Bq.  1.  Iha  quantity  k  rapraaanta  tba  diatanca 

tba  wava  baa  propagatad  In  tba  Intarval  In  taraa  of  tba  fractional  diatanca  to  tba 

dlacontinulty  diatanca  L  for  an  Initially  alnuaoldal  wava.  Introduction  of  tba 

paraaatar  0  allowa  tba  eoafflelanta  b  ,  to  bo  calculatad  fro« 

n,  i 

b^  ^  ain  tod  +  nb  f(0)]|[l  +  k  f'(0)]  d0  (3) 


Tba  dlacontinulty  diatanca  L  rapraaanta  tba  diatanca  which  Initially  alnuaoldal 
wavafon  travala  in  a  diaaipatlonlaaa  nadlun  bofora  bacowing  wultlpla  valuad.  At 
thia  diatanca,  tba  wavafront  baa  an  Infinlta  alopa,  hotfovar,  tba  paak  baa  aovad 
only  tba  diatanca  of  ona  radian  with  raapaet  to  tba  aaro  polnta  of  tba  wava.  Iha 
diatanca  L  la  givan  by 

L  -  c^|2((B/2A)  +  1]  Uj(0)*f|"’'  (U) 

wharo  f  la  tba  fraquancy  of  tba  wava  and  U.(0)  la  tba  ■axlwuw  particle  velocity 
of  the  Initial  alnuaoldal  wava.  For  thia  nodal,  thia  diatanca  L  la  aufflcient  to 
daacrlba  tba  dlatortlon  taachanlaaa  and,  conaaquantly,  la  an  la^ortant  paraaatar  in 
deacrlblng  tba  raaulta. 

The  abaorptlon  coafficiant  a  of  a  wava  of  Inflnltaaiaal  aaplltuda  with  tba 
frequency  of  the  fundaaental  cowpSnant  conplately  daacribaa  the  abaorptlon  aachanlaa; 
banco,  it  la  another  l^ortant  paraaatar  for  thia  nodal.  Each  haraonlc  la  abaorbad 
eaaantlally  exponentially  for  aach  Intarval  In  aach  Iteration.  The  property  of  heat- 
conducting,  vlacoua  aedla  having  an  abaorptlon  coefficient  dependant  on  the  aquara 
of  the  frequency  la  alao  taken  into  account. 

Ualng  the  nodal  prevloualy  daacrlbed, 
calculatlona  wore  nada  ualng  a  high  apaad 
electronic  conputar.  To  accurately  retain 
the  ahapa  of  the  wavafom,  aa  nany  aa 
alxtaen  hamonlca  ware  conputed.  Iha  >—  ■' 

value  of  k,  the  length  of  aach  Itara-  U.  (0) 

tlon  waa  taken  to  ba  0.0$,  one  twantlath 
of  the  dlacontinulty  diatanca.  Iha  raaulta  .5 
are  preaanted  In  taraa  of  the  dlaanalonlaaa 
product  at  L  which  adequately  daacribaa  both 
dlatortloS  and  abaorptlon  aachanlaaa.  Ihla 
paraaatar  a  t  la  related  to  the  invaraa  of 
the  Eaynoldi  nuabar. 

Ihla  aodal  la  not  expected  to  ba  valid  for 
axtreaely  Intanaa  aound  wavaa  which  aay 
produce  periodic  ahock  fronta  but  rather  for 
aound  wavaa  of  andarata  aaplitudaa  travailing 
In  a  fairly  abaorblng  aadlua.  Iha  raaulta 
praaentad  are  within  tbaaa  llailtatlona. 

Figure  2. 
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Flguro  3. 


Flguro  4. 


Flguro  2  ohowo  typical  raaulta  froai  calculatlona  aada  froa  thla  aodal.  Iha 
aaplltuda  of  tha  particla  valoclty  U  (K)  of  tha  harmonica  ara  nomallaad  to  tha 
aaplltudo  0.(0)  of  tha  Initial  ainuaBldal  wave.  Tha  value  of  a  L  aquala  0.10  for 
thaaa  curvai.  Tha  abaclaaa,  capital  K,  maaauraa  tha  dlatanca  tna  wava  travala  In 
unlta  of  tha  dlacontlnulty  dlatanca.  It  la  to  ba  particularly  noted  that  thaaa 
calculatlona  extend  to  10  tlaea  tha  dlacontlnulty  dlatinca.  For  thaaa  calculatlona, 
the  aaplltudaa  of  tha  higher  harmonica  have  their  peak  In  tha  region  where  K  la 
between  1  and  2. 

Figure  3  ahowa  the  effect  of  the  nonlinear  propagation  on  the  dlaalpatlon  of 
the  aound  Intenalty  calculated  from  the  prevloua  curvea.  Aa  expected,  there  la  a 
large  difference  between  tha  exponential  decay  (daahed  line)  and  tha  decay  predicted 
by  thla  model  (aolid  line).  The  predicted  greater  abaorptlon  la,  of  courae,  cauaed 
by  the  higher  rate  of  abaorptlon  of  the  genaratad  harmonica.  To  deacrlbe  thla 
hlghar  rate  of  abaorptlon,  one  may  define  a  finlta  amplitude  abaorptlon  coefficient 

^*1  “  21  dx 

Thla  coefficient  deacrlbea  tha  rate  at  which  the  energy  of  the  total  wave  la  being 
dlaalpated  at  any  given  dlatanca.  One  may  aimllarly  define  a  coefficient  for  the 
energy  carried  by  the  fundamental  component  aa 

(6) 


the 


1.0 


The  abaorptlon  coefficient  at.  neaaurea  the  rate  of  loaa  of  energy  fr( 
fundaawntal  component  alone  while ^C(.  neaaurea  that  from  the  total  wave. 

Figure  h  ahowa  theae  coefflclenta  (nomallied 
to  a  )  for  the  parameter  a  L  •  0. 10.  It  la  .Impor* 
tant^to  note  that  the  abaorptlon  coefficient  a.  for 
the  fundamental  frequency  dlffara  from  the  total 
abaorptlon  coefficient  oL  for  moat  valuea  of  K. 

It  can  be  concluded  that  the  total  abaorptlon 
coefficient  cen  not  be  alwqra  obteined  by  meaaure- 
menta  of  the  fundamental  componant  only.  The  differ* 
encea  between  QL  and  a(_  can  be  very  large  at  amall 
dlatancea  whara^the  energy  carried  by  the  funda¬ 
mental  la  balng  loat  by  both  dlatortlon  and 
abaorptlon,  whlla  tha  enargy  of  the  total  wave 
la  being  loat  only  by  abaorptlon. 

Aa  the  raaulta  depend  on  the  parameter  a.L 
la  Intereatlng  to  Inveatlgate  the  dlaalpatlon''o: 
the  Intenalty  aa  a  function  of  thla  parameter. 

Figure  5  ahowa  a  aami-logarlthmlc  plot  of  the 
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Intanalty  for  differant  valuaa  of  a  h.  For  larga  valuas  of  a  t(>0.5);  whara  tha 
absorption  la  larga  or  uhara  the  nonllnaar  affacta  are  saMll,'*tha  decay  is  nearly 
exponential  as  tha  loner  curve  in  Figure  3  is  alnost  a  straight  line.  In  these 
cases,  tha  wave  is  absorbed  so  fast  that  tha  nonllnaar  effects  do  not  becotw 
appreciable.  Hence,  there  is  little  or  no  excess  absorption.  For  saallor  values  of 
Ot  L,  the  wave  beconea  appreciably  distorted  and  the  excess  absorption  baeones 
noticeable.  Tha  uppar-nost  curve  of  Figure  ?  shows  that  tha  axcass  absorption 
becoMS  proninent  after  the  wave  has  propagated  to  the  vicinity  of  tha  discontinuity 
distance. 

The  differential  absorption  coefficient  0(-  gives  a  nore  detailed  view  of  the 
rata  of  absorption.  Figure  6  shows  the  ebsorption  coefficients  corresponding  to 
the  curves  in  Figure  5.  For  the  snaller  values  of  a  L,  the  absorption  coefficients 
have  their  suuclstf  in  tha  vicinity  of  K  >  1.6  L.  This  corresponds  to  the  situation 
that  the  peak  of  the  initial  wave  s»ve  x/S  with,  respect  to  the  sero  points  of  the 
wave.  Physically,  the  wovenent  of  this  peak  is  United  to  this  range  as  the  particles 
cannot  wov  through  each  other. 

In  Figure  ^  tha  nsxlnun  values  obtained  by  the  absorntion  coefficient  are 
plotted  at  a  function  of  a  L.  For  values  of  large  at  L(>1),  the  finite  anplitude 
distortion  is  completely  negligible.  For  small  values  of  a  L,  i.e.,  large  non¬ 
linear  affects,  the  msximum  absorption  bscosws  inversely  proportional  to  the  para¬ 
meter  at.  In  other  words,  at  large  sound  pressures,  the  msximimi  absorption 
coefficient  is  directly  proportional  to  the  Initial  sound  pressure.  The  nature  of 
the  curve  in  Figure  7  agrees  somewhat  with  the  theoretical  results  as  outlined  in 
the  review  paper  of  Zarsmbo  and  Krasil'nikov  (2). 

Although  this  model  approach  is  an  over  simplification  of  the  problem  of 
flnite-asiplltude  distortion,  the  results  obtained  from  it  are  beneficial  in  under¬ 
standing  soma  of  the  properties  of  finite-amplitude  distortion.  At  present  there 
is  not  a  sufficient  amount  of  experimental  data  to  check  tha  range  of  validity  of 
Che  results  obtained  using  this  model.  (This  work  was  supported  by  the  Office  of 
Naval  Kasaarch,  U.8.  Navy  and  Che  U.S.  Arny  lesearch  Office  (Ihirham)). 
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Optical  Effects  of  Ultrasonic  Waves  Producing  Phase  and  Amplitude  Modulation 
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A  theory  is  developed  for  the  diffraction  of  light  by  ultrasonic  waves  of  sufficiently  high  frequency,  large 
amplitude,  and/or  large  beamwidth  that  the  emerging  light  wavefront  is  significantly  amplitude  modulated 
in  addition  to  the  phsM  modulation  considered  in  the  Raman-Nath  approach.  The  ultrasonic  beam  is  con¬ 
sidered  to  consist  of  N  adjacent  sections  and  the  final  diffraction  spectrum  to  result  from  N  successive 
diffractions.  The  diffraction  orders  emerging  from  a  given  section  are  considered  to  be  sources  for  further 
diffraction  by  the  next  section.  Only  phase  noodulation  of  the  separate  plane  waves  (diffraction  orders)  is 
considered  within  a  section.  Refraction  of  light  is  not  considered  as  such;  it  is  characterised  by  successive 
redistribution  of  light  in  the  diffraction  orders.  Numerical  results  are  obtained  by  iterations  using  an  elec¬ 
tronic  computer.  Ihese  results  are  compared  with  measurements. 


INTRODUCTION 

The  difiraction  of  light  by  ultrasonic  waves  has 
been  the  subject  of  numerous  experimental  and 
theoretical  investigations.  Raman  and  Nath*  developed 
a  theory  for  conditions  under  which  the  ultrasonic  fre¬ 
quency,  amplitude,  and/or  beamwidth  are  sufficiently 
small  that  the  r^racfion  of  light  within  the  ultrasonic 
beam  may  be  neglected.  Under  such  conditions  the 
ultrasonic  beam  may  be  considered  to  act  as  a  pure 
phase  grating  producing  only  changes  in  the  relative 
phase  of  the  initially  plane  incident  light  wave  and  no 
intensity  changes,  llie  Raman-Nath  approach  has 
proven  useful  over  a  limited  range.  Herein,  ^e  Raman- 
Nath  theory  and  other  theories  based  on  the  phase 

*  Present  address:  Bell  Telephone  I-sboratories,  Inc.,  Murray 
HiU,  New  Jersey. 

•  C.  V.  Raman  and  N.  S.  Nath,  Proc.  Indian  Acad.  Sd.  A2, 406- 
412  (1935);  A3,  75-84  (1936). 


grating  concept  are  called  the  Raman-Nath  approach 
for  sake  of  brevity. 

When  the  ultrasonic  frequency,  amplitude,  and/or 
beamwidth  become  large,  refraction  causes  significant 
amplitude  modulation  along  the  emerging  light  wave- 
front.  Extermann  and  Wannier,'  Wagner,'  Van  Cittert,' 
and  Mertens*  have  obtained  solutions  for  such  condi¬ 
tions.  Their  results  contain  varying  degrees  of  approxi¬ 
mation  and  complerdty. 

In  this  paper,  a  solution  to  the  problem  of  difiraction 
of  light  by  sinusoidal,  plane,  progressive,  ultrasonic 
waves  is  presented.  The  ultrasonic  beam  is  considered 
to  consist  of  N  adjacent  ultrasonic  beams.  For  N 

'  R.  Extermann  and  G.  Wannier,  Helv.  Phys.  Acta.  9,  520-532 
(1936). 

•  E.  H.  Wagner,  Z.  Phyaik  141, 604-621  (1955). 

*  P.  H.  Van  attert,  Phydca  4,  (1937\ 

'  R.  Mertens,  Meoedef.  Koninkl.  Viaam.  Acad.  Wetenschap. 
Belg.  12, 1-37  (1950). 
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Fig.  1.  Schematic  dUjiram  showing  a  diffraction  order  with 
amplitude  ♦,<•“‘>(0)  from  the  (a— 1)  st  section  ii^dent  at  *— 0 
on  the  nth  section  with  angle  of  incidence  d,;  sound  propagation 
in  the  +x  direction,  incident  light  in  the  +s-direction. 

sufficiently  large  only  phase  modulation  of  the  light 
wavefronts  need  be  considered  within  each  of  the  N 
intervals.  The  final  diffraction  spectrum  is  considered 
to  result  from  N  successive  diffractions. 

THEORY 

Consider,  as  indicated  in  Fig.  1,  a  section  within  the 
ultrasonic  beam  having  thickness  I  assumed  to  be 
sufficiently  small  that  refraction  of  light  within  the 
distance  I  may  be  neglected.  Let  this  be  the  Mth  of 
many  identical  sections  making  up  the  whole  ultrasonic 
beam.  Diffraction  orders  emerging  from  the  (»— l)st 
section  may  be  considered  sources  for  further  diffraction 
by  the  »th  section.  Figure  1  schematically  represents 
the  following  situations :  A  diffraction  order  with  ampli* 
tude  front  the  («— l)st  section  is  mcident  at 

x=0  on  the  nth  section  with  angle  of  incidence  0,. 
Neglecting  refraction,  the  incident  light  undergoes  a 
change  in  phase  on  progressing  through  the  nth  section. 
This  change  is  expressed  by* 

[-2ir»  /•"««*»  T 

—  J  (1) 

where  the  integral  in  Eq.  (1)  represents  the  optical 
path  through  the  nth  section,  X  is  the  wavelength  of 
light  in  vacuo,  and  fi(r)  expresses  the  refractive  index 
along  the  light  path.  Using 

lt(.s)=ito+U  sin2TC(j  sinff,-*')A*]i  (2) 

where  X*  is  the  wavelength  of  sound  and  the  Raman- 
Nath  parameter  Vi = (3) 
Eq.  (1)  becomes 


'In  this  development  the  time  dependence  of  the  light  ud 
sound  are  not  included,  with  the  result  that  the  Dialer  shifts  of 
the  various  diffraction  orders  ate  not  shown.  The  small  Doppler 
shift  is  neglected  in  writing  the  coefficient  of  the  integral  in  Eq.  (1). 


Using  the  identity 
exp(— ta  sini) 

■He 

“exp[ta8in(— ft)]=  £  /,(o)  exp(— 156),  (S) 


Eq.  (4)  may  be  expressed  in  the  form 


In  Eq.  (6)  a  term  for  a  particular  q  represents  light 
incident  on  the  »th  section  at  an  angle  6,  and  incident 
on  the  (n+l)st  section  at  an  angle  where 

sinff«=— kX/(poX*).  (7) 

Consider  one  such  component  (particular  p  and  9)  from 
Eq.  (6)  denoted  by  0,, ,'"'(**)•  This  relates  a  particular 
emerging  component  at  x'  to  the  one,  namely  (0), 
incident  on  the  «th  section  at  *=0.  In  order  ffiat  each 
successive  section  may  be  treated  in  mathematically 
identical  manners,  roust  be  transformed  to 

in  order  that  the  point  of  incidence  on  the 
(»i+l)st  section  will  also  be  at  *=0. 

The  transformation  from  »'  to  involves  an 
ultrasonic  phase  shift  and  an  optical  phase  shift  as  indi¬ 
cated  in  Fig.  2.  The  effect  of  an  ultrasonic  phase  shift 
is  expressed  by  a  factor  exp(— where  6*  is  the 
ultrasonic  phase  difference.  From  Fig.  2  it  is  seen  that 

exp(— t^i*)  =  exp(2irt^/ tanff,/X*).  (8) 

Similarly,  the  effect  of  the  optical  phase  shift  is  ex¬ 
pressed  by  a  factor 

exp(— t8)=exp(2irt#»o/tanff,  sinff,/X) 

=exp(— 2irt^f  tanff,/X*), 

where  i  is  the  optical  phase  difference  and  Eq.  (7)  has 
been  used.  The  transformation  then  consists  of  multi- 


Fio.  2.  Schematic  diagram  illustrating  the  optical  and  ultra¬ 
sonic  phase  differences  involved  in  transforming  from 
to  X— O. 
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plying  by  the  right-hand  members  of  Eqs. 

(8)  and  (9).  Noting  that  these  two  factors  are  simply 
complex  conjugates,  and  that  their  product  is  unity, 
then 

(10) 

Having  transformed  the  separate  components  corre¬ 
sponding  to  various  p  and  9  to  a  common  point  x^O, 
the  components  propagating  in  a  given  direction  may 
be  combined  to  obtain  the  amplitudes  of  the  diffraction 
orders  which  act  as  new  sources  for  further  diffraction 
in  the  (»-l-l)st  section.  The  terms  to  be  combined  to 
obtain  the  amplitude  of  the  rth  order  are  those  for 
which  p+q=f,  i.e.. 


+•0 

'rl  tanff,\T 

•'>  sm( 

- -)\ 

L/t  smBp  > 

.  X*  /J 

'-2TtMo/l  rirt(r 

-  P)l 

(11) 

Xexp 

exp 

.  X  cosB,  J  L 

X.  ]■ 

where  x'  has  been  expressed  by 

x'=l  tanflp. 

(12) 

Equation  (11)  expresses  an  iterative  procedure  for 
calculating  the  final  diffraction  order  amplitudes 
resulting  from  N  successive  diffractions. 

Equation  (11)  may  be  put  into  an  approximate  form 
which  is  more  suitable  for  numerical  evaluation.  Using 
the  approximations 

tanfl,2isinSt  =  —  k\/  (moX*)  ,  (13) 

l/cosfit=sec5*£^l-|-j9t’,  where  fft^sinfl*,  (14) 


Fio.  4.  Intensity  modulation  of  the  light  wavefront  emer|^  from 
the  ultrasonic  beam  for  vr**4  and  Qr- 1.5,  using  Qt'. 


In  either  Eq.  (11)  or  Eq.  (16) 


W“>  =  1, 

and 


(17) 

(18) 


for  unit  incident-light  amplitude  and  normal  light  inci¬ 
dence.  The  light  intensities  in  the  final  spectrum  are 
obtained  from 


(19) 


It  has  been  assumed  that  the  refraction  of  diffraction 
orders  within  each  section  may  be  neglected  and  the 
Raman-Nath  approach  appli^  therein.  Limitations 
on  Q  and  v  for  validity  of  the  Raman-Nath  approach 
have  been  given  in  the  form* 


and  introducing  the  parameter' 


ei«2. 


(20) 


g,=  (2TX/)/(MoX«), 

Eq.  (11)  becomes 

p- 


(15) 


\  iPQi  / 


Xexp(-|*^i).  (16) 


Fio.  3.  Light  intensities  predicted  by  the  present  theory 
(vertical  bars)  for  Qr~1.5,  using  Qt'  in  Eq.  (16)  and  by  the 
Raman-Nath  theory  (circles}. 


’  The  parameter  Q  ^>pears  in  some  form  in  various  other 
theories  tor  the  diffraction  of  light  by  ultrasonic  waves. 


The  degree  of  inequality  required  in  (20)  depends  on  the 
accuracy  required.  It  is  proposed  that  the  number  of 
sections  chosen  for  application  of  the  present  theory  be 
suflSciently  large  that  (0p)i  be  less  than  two  by  a 
factor  of  one  or  two  orders  of  magnitude.  The  number 
of  sections  N  may  be  increased  as  required  to  obtain 
theoretical  results  for  large  Qt  and  vr  pertaining  to  the 
total  ultrasonic  beam. 

Straightforward  application  of  the  present  theory 
would  require  that 

I^L/N,  (21) 

where  L  is  the  total  width  of  the  ultrasonic  beam,  be 
used  in  Eqs.  (3)  and  (15)  which  define  vi  and  Qt  which 
appear  in  Eq.  (16).  However,  an  alternative  approach 
is  proposed.  Note,  in  Eq.  (18),  that  in  the  first  section 
the  result  is  independent  of  Q.  Refraction  effects,  as 
characterized  by  spreading  of  light  into  the  various 
difraction  orders  as  it  emerges  from  each  succes.sive 
section,  are  in  no  way  considered  in  the  first  section. 
If  the  ultrasonic  beam  width  L  were  divided  into  twice 
as  many  sections,  refraction  effects  would  be  ignored 
in  a  first  section  of  half  the  width.  In  order  to  approxi- 
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Fig.  S.  Theoretical  (line)  and  ezperimentai  (circles)  zeroth- 
order  light  intensity  vs  Raman-Nath  parameter  for  ^"0.62, 
using  Qt'. 


mately  account  for  refraction  eflfects  in  a  first  interval 
of  finite  width,  it  is  proposed  that 

0/=Cr/(iV-l)  (22) 

be  used  in  Eq.  (16)  to  replace  Qt.  The  Raman-Nath 

parameter  tt  should  remain  as 

v,=tT/N.  (23) 

Use  of  Eq.  (22)  is  equivalent  to  using 

Qt'=QtN/(N-1)  (24) 

as  an  “effective  Q”  to  compensate  for  use  of  finite 
sections  of  the  ultrasonic  beam.  Obviously,  as  the 
number  of  sections  N  becomes  very  large,  the  difference 
between  Qt'  and  Qt  vanishes. 

mrMERICAL  RESTTLTS  AND  COMPARISON  WITH 
EXPERIMENTAL  DATA 

Numerical  calculations  were  made  from  Eq.  (16) 
using  the  Michigan  Sute  Integral  Computer  (MISTIC). 
The  choice  of  0<»r<8,  0r<1.5,  and  iV=16  gives 
(Qi))j<0.05.  The  value  A^=  16  was  also  chosen  for  con¬ 
venience  in  a  binary  computer.  Calculations  were  made 
for  0.25  intervals  of  vt  and  for  p  (and  therefore  also  r) 


Fig.  6.  Theoietlcal  6r»t-order  light  intensity  vs  Rsman-Nnth 
parameter  for  Qr**0.62,  using  Qt'. 


ranging  from  —4  to  +4.  For  a  given  range  of  vt  the 
range  of  p  should  be  taken  over  values  for  which  signifi¬ 
cant  light  occurs  in  the  ^th  orders.  The  final  results  are 
probably  less  reliable  as  |r|  approaches  the  maximum 
1^1  used. 

The  differences  between  predictions  of  the  present 
theory  and  the  Raman-Nath  theory  are  illustrated  in 
Fig.  3.  The  predicted  light  intensities  in  the  zeroth 
through  fourth  diffraction  orders  are  shown  for  Qt=  l-S, 
using  Qt  in  Eq.  (16),  for  »r=2,  4,  and  6. 

Using  the  calculated  amplitudes  and  phases  of  the 
diffraction  orders,  the  calculated  intensity  modulation  of 
the  light  wavefront  emerging  from  the  ultrasonic  beam 
for  the  Hr =4  case  shown  in  Fig.  3  is  shown  in  Fig.  4. 
Note  that  the  intensity  modulation  is  quite  significant. 
This  intensity  modulation  explains  the  failure  of  the 
Raman-Nath  theory  for  this  case  where  (Q»)r=6. 
Taking  the  refraction  viewpoint,  concentration  of  light 
about  Jt  indicates  light  refracted  toward  this  region 
where  the  density  of  the  medium  and  hence  the  refrac¬ 
tive  index  is  greatest. 


Fig.  7.  Theoretical  and  ezperimentai  zeroth.order  light  intenrity 
VB  Raman-Nath  parameter  for  Or“0.93,  using  Qt'. 


Figures  5  through  10  show  calculated  zeroth-  and  first- 
diffraction-order  light  intensities  for  0r=O.62,  0.93, 
and  1.24,  using  Qt  in  Eq.  (16).  The  predictions  of  the 
Raman-Nath  theory  are  also  shown  in  these  figures, 
where  clarity  permits.  Experimental  results  for  the 
indicated  values  of  Qt,  obtained  by  Klein,*  are  also 
shown  in  Figs.  5  and  7  through  10.  Klein’s  measure¬ 
ments  were  made  at  5.23  Me  in  water.  Qt  was  varied  by 
varying  the  ultrasonic  beam  width  (2.0, 3.0,  and  4.0  cm 
for  the  data  shown). 

Klein  found  good  agreement  between  Mertens’* 
theory  and  his  experimental  values  of  light  intensity 
in  the  zeroth  and  first  diffraction  order,*  for  a  limited 
range  of  Qt  and  vt-  The  present  theory  agrees  with 
the  Mertens  theory  where  the  Mertens  theory  is  in 
reasonable  agreement  with  measurements.  However, 
the  present  theory  appears  to  be  valid  over  a  wider 
range  of  Qt  and  Vt. 

•  W.  R.  Xtein,  M.  S.  thesis,  Michigan  State  University  (1962). 
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For  comparison,  the  zeroth-order  light  intensities 
predicted  by  the  present  theory  (using  both  Qr>»1.5 
andOr'<“l.S(Ar/JV— 1)  for  calculations),  by  theMertens’ 
theory  for  Qr=1.5,  and  by  the  Raman-Nath  theory 
are  shown  in  Fig.  11.  If  the  present  theory  is  valid,  then 
the  Mertens  theory  is  a  good  approximation,  in  this 
case,  for  vr<4.0.  The  difference  between  the  curves  in 
Fig.  11  for  the  present  theory  using  Qt  and  Qt  is  not 
great,  but  the  difference  is  generally  in  the  correct 
sense  to  give  better  agreement  with  experimental  values 
by  using  Qt'. 

Calculations  have  been  made  for  Qr=0.31  and  com¬ 
pared  with  the  predictions  of  the  Raman-Nath  theory. 
The  predicted  light  intensities  in  the  diffraction  orders 
agree  to  within  0.01  up  to  the  values  of  vr  indicated  in 
Table  I.  This  agreement  over  a  fairly  large  range  of  t>r 
demonstrates  that  the  method  of  calculation  used  in 
the  present  theory  gives  results  which  approach  the 
Raman-Nath  results  for  a  small  value  of  Qt. 

It  should  be  pointed  out,  however,  that  only  agree¬ 
ment  in  light  itUensily  has  been  obtained  for  ^t=0.31. 
The  relative  optical  phases  should  also  be  considered. 
The  intensity  modulation  of  the  emerging  light  wave- 
front  shown  in  Fig.  4  for  Qt=  1.5  results  not  only  from 
light  amplitudes  which  differ  from  those  predicted  by 


Fio.  8.  Theoretical  and  experimental  first-order  light  intenrity  vs 
Raman-Nath  parameter  for  Qr-0.93,  using  Qt', 


Fio.  9.  Theoretical  and  experimental  zeroth-order  light  intensity 
vs  Raman-Nath  parameter  for  Qt“  1.24,  using  Qt’. 


Fig.  10.  Theoretical  and  experimental  first-order  light  intensity 
vs  Raman-Nath  parameter  for  Qr’*1.24,  using  Qt'. 


the  Raman-Nath  theory,  but  also  from  different  rela¬ 
tive  phases  of  the  various  diffraction  orders.  The  relative 
phases  in  the  Raman-Nath  theory  are  either  0  or  r  rad, 
according  to  whether  the  Bessel  function  is  positive  or 
negative,  respectively.  The  calculations  for  0r=O.31, 
though  the  light  intensities  agree  with  the  Raman-Nath 
theory,  show  marked  deviations  in  relative  phase  from 
0  or  T  rad.  This  indicates  that  the  Raman-Nath  theory 
may  appear  satisfactory  from  intensity  measurements 
but  not  be  accurate  where  the  relative  phases  are  con¬ 
cerned.  The  relative  phases  become  important,  as 
pointed  out  in  a  previous  paper,'  in  a  situation  in¬ 
volving  successive  diffraction  of  light  by  two  separate 
ultrasonic  beams. 


Table  I.  Maximum  values  of  ir  for  which  the  predictions  of 
the  present  theory  and  the  Raman-Nath  theory  agree  to  within 
0.01  light  intensity.  Qr~0.31. 


Order 

0 

±1 

±2 

±3 

±4 

VT 

7.0 

>8 

7.25 

6.0 

4.0 

'  L.  E.  Hargrove,  E.  A.  Hiedemarm,  and  R.  Mertens,  Z.  Physik 
167,  326-336  (1962). 
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Fig.  12.  Intensity  modulation  of  the  light  wavefront  emerging  from 
the  ultrasonic  beam  for  Fi— 4  and  Qr— 0.31,  using  Qt'. 


The  agreement  in  intensity  predictions  may  also  be 
considered  to  demonstrate  that  the  Raman-Nath 
approach  predicts  the  correct  intensities  even  when 
tWe  is  some  degree  of  light  intensity  modulation  along 
the  emerging  light  wavefront.  Figure  12  shows  the  in¬ 
tensity  modulation  calculated  for  Qr=0.31  and  »t=4.0. 
It  has  been  shown  in  Table  I  that  the  present  theory  and 
the  Raman-Nath  theory  give  essentially  the  same  light 
intensities  through  fourth  orders  up  to  this  value  of  vt. 
Therefore,  the  intensity  modulation  shown  in  Fig.  12 
is  an  estimate  of  the  upper  limit  of  intensity  modulation 
for  which  the  Raman-Nath  theory  gives  the  correct 
light  intensities  in  the  first  four  diffraction  orders. 
Figures  4  and  12  are  for  the  same  value  of  vr  but 


different  Qt,  and  may  be  compared  to  illustrate  the 
difference  in  intensity  modulation  for  ultrasonic  beams 
with  the  same  “diffraction  strength”  vt  but  different 
Qt-  Finally,  it  is  remarked  that  the  light  intensity 
modulation  for  the  smallest  (Qt)r  calculated  (Qr=0.31 
and  iir=0.25,  giving  (Q»)r~0.08)  deviates  from  unity 
by  at  most  approximately  0.04. 

DISCUSSION 

The  predictions  of  the  present  theory  are  in  good 
agreement  with  Klein’s  experimental  data.  The  agree¬ 
ment  is  somewhat  better,  for  the  larger  values  of  Qt 
and  hr,  than  that  found  by  Klein  using  Mertens’  theory. 

Diffraction  of  light  by  high-frequency,  intense,  ultra¬ 
sonic  waves  is  of  particular  interest  for  determining 
the  waveform  of  distorted  finite-amplitude  waves  in 
liquids.  While  the  present  results  are  restricted  to 
sinusoidal  waves,  work  is  in  progress  to  extend  the 
theory  to  include  arbitrary  ultrasonic  waveform, 
especially  waveforms  of  the  finite-amplitude  type. 
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Varioua  anparlBMntal  taehniquas  for  tha  production  of  aurfSca  wavaa  on  a  fraa 
aurfaea  of  a  solid  hava  baan  dascribad  in  tba  past  (1).  Quits  fraguantly  ona  uses 
tba  nathod  of  convarting  a  longitudinal  uanra  in  a  aolid  to  a  aurfaea  wave  on  another 
solid  in  contact  with  tha  first,  lha  prasant  paper  ia  eoneamad  with  this  type  of 
■ode  eonreraion  because  the  energy  transfer  frcai  a  longitudinal  to  a  surface  wave 
seans  to  offer  poaalbillties  for  decreasing  the  intensity  of  the  reflection  of  the 
longitudinal  wave  in  the  aolid  through  which  it  la  propagated.  Sena  eaperlnental 
results  for  a  specific  ease  are  given  here. 

The  production  of  a  surface  wave  is  usually  aceo^lished  by  the  fanlllar  wedge 
nethod  ahoun  in  Figure  1.  The  anplituda  of  the  surface  wave  reaches  a  aMuriaun  (2) 

whan  tha  angle  of  incidance  a  is  chosen  that  the  ratio 
Vt,/Ve  m  aino,  idiere  Vj,  is  the  velocity  of  the  longitu* 
dlnal  wave  in  solid  I,  and  Vg  ia  the  velocity  of  the 
surface  wave  on  solid  II.  An  arrangaawnt  of  this  type 
has  the  disadvantage  that  ona  cannot  easily  detemlna 
bow  wuch  energy  contained  in  the  longitudinal  wave  is 
transferred  to  the  surface  wave.  A  slightly  different 
wedge  was  therefore  used,  shown  ia  Figure  2.  An  ultra* 
aonic  pulse  is  awlttad  by  transducer  B.  If  the  botton 
surface  of  the  wedge  is  in  contact  with  air  the  pulse ' 
will  be  totally  reflected  and  will  ba  picked  up  by  the 
recalving  transducer  R.  Sie  aavlltude  of  the  received  pulse  is  taken  as  reference 
for  the  subsequent  eaparlnent  in  which  the  wedge  is  placsd  on  the  anooth  surface  of 
solid  II.  In  the  eaparisNats  described  here  solid  II  is  alualnuai  while  solid  I  is 

Flaalglas.  Ao  velocities  Vj^  and  Vg  were  deteminad 
axperinantally  by  a  conventional  pulse  awthod.  Flee* 
ing  the  Flexiglaa  wedge  on  the  aluninun  surface  allows 
the  nods  convarsion  froa  a  longitudinal  wave  in  tha 
wedge  to  a  surface  wave  on  the  aluninun.  Ae  wedge  is 
cut  in  such  a  way  that  the  angle  of  incidence  a  satis* 
fiss  ths  abovs  condition  for  optlann  aurfaea  wave  pro¬ 
duction.  With  the  wedge  on  solid  II  there  is  no  nore 
total  reflection  of  the  incident  pulse,  and  the  anpli- 
tude  of  the  signal  received  by  tranaducer  K  should  de¬ 
crease  by  an  anount  detacnined  by  the  energy  converted 
to  the  surface  wave.  However,  it  is  found  that  using  tha  Flexiglaa  wedge  the  anpli- 
tuda  of  the  received  pulse  increases  idtan  the  wedge  ia  placed  on  the  aluninun.  Ais 
nay  ba  explained  by  considering  the  loss  nechanisna.  If  the  wedge  is  not  in  contact 
with  solid  II  every  part  of  the  sound  signal  anittod  by  tranaducer  B  travels  the  sane 
distance  in  the  wedge  before  it  reaches  transducer  E.  Ae  absorption  in  tha  Plexi¬ 
glas  is  quite  high.  If  the  wedge  is  than  placed  on  the  aluninun  a  surface  wave  is 
generated  between  tha  two  solids.  As  this  wavs  propagates  along  the  surfacea  it  nay 
radiate  energy  back  into  the  wedge.  If  the  absorption  of  this  wave  should  be  lass 
than  the  absorption  of  the  longitudinal  wave  in  the  wedge  a  signal  of  higher  aapli- 
tuds  would  be  received  at  R  than  if  tha  signal  had  traveled  entirely  in  tha  Plexiglas. 
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Om  would,  thwraforu.  Ilk*  to  ■InlmlM  thla  radiation  back  Into  aolld  I.  Idaally, 
on*  would  Ilk*  to  aoparat*  th*  aurfae*  wav*  frow  th*  bottoa  of  th*  w«dg*  aa  aoon  aa 
It  la  ganaratad.  But  thlt’wav*  la  ganaratad  ovar  a  flnlta  araa  of  contact  batwaan 
th*  two  aolida,  and  th*  ala*  of  thla  araa  la  datarwlnad  by  th*  ala*  of  tha  awlttlng 
tranaducar  and  th*  angl*  of  Ineldaneat  on*  tharafor*  ancountara  eartaln  aaparlaantal 
l^tatlona  In  aaparatlng  tha  aurfae*  wav*  frow  th*  bottca  of  th*  wadg*.  Sine*  th* 
cantral  portion  of  tha  ultraaonle  algnal  awlttad  frow  tranaducar  E  contain*  wor*  anar- 

gy  than  th*  outalda  of  th*  baaa  on*  alght  rawov*  th* 
oorraapondlng  woat  Intana*  portion  of  thp  aurfae*  wav* 
by  rounding  off  on*  adg*  of  aolld  II  and  placing  It  In 
ralatlon  to  th*  wadg*  aa  Indleatad  In  Flgura  3.  Ih* 
radlua  of  curvatur*  of  tha  roundad  adg*  la  big  anough 
ao  that  tha  aurfae*  wava  will  traval  around  It  (3). 
Adjuatlng  th*  ralatlv*  poaltlona  of  th*  wadg*  and  aolld 
II  In  thla  wannar  on*  traaafar*  tha  aoat  Intana*  part 
of  th*  longitudinal  wav*  In  tha  wadg*  to  a  aurfae*  wave 
which  la  than  conducted  away  frea  th*  wedge,  thua  re¬ 
ducing  radiation  back  toward  th*  receiving  tranaducar. 

In  ordar  to  dacraaa*  further  th*  are*  of  contact  In  which  aurfae*  wava*  are  gana¬ 
ratad  on*  can  dacraaa*  th*  dlaanalona  of  tha  aadttlng  tranaducar.  Doing  thla  on*  nar¬ 
row*  tha  width  of  tha  longitudinal  wav*  Inplnglng  on  aolld  II  and  on*  can  reduce  th* 

received  algnal  by  replacing  th*  alunlnun  block  by  a 
aolld  alunlnun  cylinder  aa  ahown  In  Flgura  k.  Thla  ra- 
dueaa  th*  araa  of  contact  batwaan  th*  two  aolida  and 
aaauraa  a  nara  efficient  aaparatlon  of  th*  aurfae*  wav* 
fren  th*  botten  of  th*  wadg*.  Changing  tha  ala*  of 
th*  aluadnun  cylinder  change*  th*  affective  araa  of 
contact  and  thua  th*  anplltuda  of  th*  racalvad  pula*. 

Ih*  angl*  of  Incidence  a  la  not  affactad  by  th*  ala* 
of  th*  cylinder. 

Ih*  aurfae*  wav*  traveling  around  the  cyllndar  la 
abaorbad  by  placing  aon*  aultabl*  llguld  on  th*  cylln¬ 
dar;  If  thla  la  not  dona,  the  wav*  ccnplataa  th*  elrel* 
and  raachaa  th*  wadg*  again  whar*  It  la  converted  back 
Into  a  longitudinal  wav*  idtlch  travel*  toward  E.  In 
thla  caa*  on*  haa  a  abort  dalay  lln*. 

Ih*  efficiency  of  nod*  convaralon  froai  longitudinal  to  aurfae*  wav*  for  th*  varloua 
arrangaaanta  daaerlbad  can  ba  aaan  txcm  the  raaulta  of  axparlaanta  ualng  Plaoclglaa  and 
aluwlnua.  Ualng  a  $  Me  pula*  frew  a  1.2  x  2.$  ca  tranaducar  and  placing  tha  antlr* 
wadg*  on  tha  alualnua  on*  raealvaa  a  algnal  whoa*  awplltud*  la  dtown  In  Figure  3b. 

Iba  rafarane*  aaplltud*  (wadg*  In  air)  la  ahown  In 
Figure  3a.  Ih*  aaplltud*  of  th*  algnal  haa  Ineraaaad 
by  about  30  pareant;  thua,  placing  th*  antlr*  wadg*  on 
tha  aluwlnwB  block  doaa  not  raduc*  th*  raflaetlona  back 
Into  th*  wadg*.  Figure  3c  ahowa  th*  received  aaplltud* 
of  th*  reflected  longitudinal  wav*  recorded  with  th* 
arraagawant  Indleatad  In  Figure  3.  In  thla  caa*  th*  aa- 
plltud*  haa  daeraaaad  acawidiat.  A  aor*  affaetlv*  ener¬ 
gy  traaafar  la  attalnabl*  if  th*  wadg*  la  placed  on  a 
aolld  alualnua  cylinder.  Ih*  racalvad  aaplltud*  la  In¬ 
dicated  In  Figure  6b,  coaparad  with  th*  rafarane*  aapll¬ 
tud*  In  Figure  6*  (w^a  In  air).  Ih*  dlaaatar  of  tha 
cylinder  waa  I.83  ca.  Iha  affaetlv*  araa  of  contact 
batwaan  th*  two  aolida  waa  Ineraaaad  allghtly  by  ualng 
a  aolld  alualnua  cyllndar  with  a  dlaaatar  of  6.3k  ca. 

In  thla  caa*  th*  racalvad  algnal  had  an  aaplltud*  ahown 
In  Figure  6c.  Ih*  raaulta  ahown  In  Figure  6  war*  ob- 
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tainad  by  ualag  an  naltting  traaaduear  with  dlaanalona  0. 2?  x  2.  ^  cn  nhich  alnlnlaad 
furthar  tha  affaeclva  araa  of  aitrfaeo  aawa  production.  It  la  aaon  that  tha  rafloetad 
longitudinal  wava  daexaaaaa  to  about  half  Ita  aaplltuda  In  tha  laat  case  ahown  In 
Flgura  6. 

Ihla  reduction  In  anplltuda  of  tha  rafloetad  longitudinal  nava  Indleataa  that  an 
approelabla  part  of  tha  anargy  In  tha  longitudinal  wva  can  ba  takan  out  of  aolld  I 
by  noda  convaralon  to  a  aurfaea  uava  on  aolld  II,  aapoclally  If  ona  conaldara  that  tha 
graataat  part  of  tha  eroaa-aaetlon  of  tha  aound  baaa  la  laplnglng  on  tha  Plaxlglaa-alr 
Intarfaca  and  only  a  anall  part  on  tha  Plaatlglaa-alualnua  Intarface.  Tha  raaultlng 
aurfaea  uava  can  ba  abaorbad  vary  aaally  ao  that  It  will  net  ba  rafloetad  back  Into 
aolld  I.  That  fraction  of  tha  longitudinal  wave  In  aolld  X  which  waa  not  convortod  to 

a  aurfaea  wave  at  tha  boundary  can  ba  raduead  furthar 
by  rapaatlng  the  proeaaa  aa  Indleatad  la  Figure  7.  In 
prlnelpla.  It  ahould  ba  poaalble  to  uaa  thla  tactolqua 
for  tha  aurpraaalon  of  raflaetlona  of  aa  ultraaonle 
wava  In  a  aolld.  Since  tha  abaorptlon  of  tha  raaultlng 
aurfaea  wava  praaanta  no  problwa  ona  ahould  ba  able  to 
approach  tha  production  of  prograaalva  ultraaonle  wavaa 
In  aollda  ualng  tha  nathod  outlined  bare,  provided  tha 
velocity  of  tha  longitudinal  wava  In  tha  aolld  In  quaa* 
tlon  la  Mollar  than  tha  velocity  of  tha  aurfaea  wave 
on  aolld  II.  (Ihla  work  waa  aupportad  by  tha  Office  of 
naval  Kaaaarch,  D.S.  Havy. ) 
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Reflection  and  Refraction  of  Mechanical  Waves  at  Solid-Liquid  Boundaries 

by 

Walter  G.  Mayer 


Abstract. 

The  energy  ratios  of  reflected  and  refracted  waves  to  the 
Incident  wave  at  ten  liquid-solid  boundaries  are  calculated  as  a  function 
of  the  angle  of  Incidence.  The  influence  of  the  wave  velocities  In  the 
media  and  their  density  on  the  shape  of  the  curves  Is  discussed. 


INTRODUCTION 

Knowledge  of  the  energy  ratios  between  incident,  reflected,  and 

refracted  mechanical  waves  at  Interfaces  formed  by  two  dissimilar  media 

l.o  1*. 

Is  of  great  Importance  in  ultrasonics  ,  seismology  ,  and  for  material 

56  7  8 

testing  techniques^'  .  Knott'  and  Zoepprits^  have  given  expressions 
describing  the  changes  of  amplitude  of  reflected  and  refracted  waves 

4 

ae  a  function  of  the  angle  of  Incidence.  Based  on  their  work,  Ergln 

has  calculated  these  amplitude  ratios  for  seismic  waves  incident  at  both 

■3 

sides  of  the  ocean  floor.  Mayer  and  Kelsey'  have  measured  the  velocities 
of  ultrasonic  waves  In  solids  by  observing  the  amplitudes  of  the  reflected 
waves  In  the  liquid.  These  measurements  and  the  results  given  here  allow 
one  to  draw  a  number  of  useful  conclusions  about  the  behavior  of  the  re¬ 
flected  and  refracted  waves.  However,  Ergln  considers  only  the  changes 
in  the  energy  ratios  caused  by  a  relatively  small  change  In  the  ratio  of 
the  longitudinal  to  shear  wave  velocity  (essentially  the  Poisson's  ratio) 
of  the  solid  medium.  The  present  paper  also  considers  the  Influence  of 
the  density  ratio  and  the  ratio  of  the  longitudinal  wave  velocities  of 
the  two  media. 
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RESULTS 

The  syetem  under  eoneideratlon  conelata  of  a  liquid  with  denaicy 
and  a  aolld  medium  with  denaity  p^.  The  plane  aurface  of  the  solid 
la  in  contact  with  the  liquid.  The  velocity  of  the  longitudinal  wave  in 
the  liquid  ia  The  velocity  of  the  longitudinal  wave  in  the  eolid  ie 

given  by  and  that  of  the  ehear  wave  by  The  energy  of  the  incident 

wave  in  the  liquid  la  unity,  and  the  angle  of  Incidence  Is  denoted  by  a, 

4 

measured  from  the  normal  to  the  liquid-solid  boundary.  According  to  Ergln  , 
the  energy  ratio  of  reflected  to  incident  wave  is  given  by 

(R/I)^  b  ^[cos  0  -  Aco.sa(l  -  B)]/[cos  a  +  Acosck(1  -  b)}^  (1) 

where  A  =  Vj^gPg/V^^jP^, 

B  =  ZBiny9in2y[coBy  -  (Vgg/Vj^2)cos3], 

The  angles  3  and  y  are  the  angles  of  refraction  of  the  longitudin¬ 
al  and  the  shear  wave  In  the  solid,  determined  by  Snell's  law 

Vj_j/slna  =  Vj^2/sin3  =  Vg^/slny. 

The  energy  ratio  of  refracted  longitudinal  wave  in  the  solid 
to  the  incident  wave  in  the  liquid  is  given  by 

(L/l)^  -  ^2cos27  (Acosa  cosy)?/  [cosP  +  Acosci((l  -  B)]^  .  (2) 

The  energy  ratio  of  refracted  shear  wave  to  incident  wave  is  then 

(S/I)2  -  1  -  (R/1)^  -  iUlf.  (3) 

Calculations  of  these  ratios  were  made  for  a  number  of  solids 
with  the  Incident  wave  In  water  and  oil.  The  values  of  velocities  and 
densities  used  are  shown  In  Table  1.  Poisson's  ratio  a  is  also  shown. 
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Table  X.  Veloeltiea  and  Denaltlea  Uaed  for  Calculations 


\l 

B 

P2 

a 

Water 

IU90 

1.00 

Oil 

17h0 

0.87 

Steel 

5850 

3230 

7.80 

0.281 

Brass 

hh30 

2123 

8.10 

0.351 

Copper 

4700 

2260 

8.90 

0.350 

Aluminum 

6330 

3I3D 

2.70 

0.338 

Magnesium 

5770 

3050 

1.70 

0.306 

Velocities  in  m/sec. 

Since  one  cannot  readily  see  from  eqs.  (1*3)  bow  the  Intensity 


of  the  waves  Is  Influenced  by  changes  In  the  various  parameters.  Individual 
calculations  were  made  for  all  the  possible  combinations  of  liquids  and 
solids  listed  In  Table  I.  The  results  are  plotted  In  Figs.  1*%  The 
flolld  lines  are  for  oil- sol Id  interfaces,  the  dashed  lines  for  water- 
Molld  Interfaces.  The  letters  R,  L.  and  S  denote  the  curves  for  the 
reflected  longitudinal  wave,  the  refracted  longitudinal  wave,  and  the 
refracted  shear  wave,  respectively. 

DISCUSSION 

The  general  conclusions  drawn  from  the  curves  In  Figs.  1-^  are 
valid  only  If  This  excludes  boundaries  formed  by  some 

liquids  and  plastics  where 

The  ratio  (R/X)^  becomes  unity  at  the  critical  angle  for  the 
refracted  longitudinal  wave  (sin^  >  1)  and  again  at  the  critical  angle 
for  the  shear  wave  (slny  •  1).  The  location  of  these  critical  angles 
0^  and  Og  depends  only  on  the  ratios  ^S2^\i*  reflection 

occurs  at  these  angles  of  Incidence  regardless  of  the  values  of  density 
and  velocity. 
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The  value  of  (R/l)  for  normal  incidence  depende  only  on 

and  eq.  (1)  reduces  to  the  well-known  formula  for  the  reflection 

coefficient.  The  shear  wave  velocity  Is  of  no  consequence  at  oe  b  0. 

2 

The  general  behavior  and  the  minima  of  the  (R/I)  curve  remains 
essentially  the  same  If  Instead  of  given  values  for  and  p^  one  selects 
slightly  different  values  without  changing  appreciably  the  acoustic  im¬ 
pedance  of  the  liquid  medium  (water  and  oil).  This  also  holds  for  the 
2  2 

maxima  of  the  (L/l)  and  (S/l)  curves. 

Poisson's  ratio  Is  not  a  dominating  factor  In  the  behavior  of 
the  curves.  Of  the  solids  considered  steel  and  magnesium  have  the  lowest 
Poisson's  ratios,  yet,  the  curves  for  the  other  three  substances  lie  between 
those  for  steel  and  magnesium. 

The  energy  of  the  reflected  wave  is  determined  mainly  by  its 

2 

value  at  a  B  0.  The  ratio  (R/l)  stays  almost  constant  until  the  critical 
angle  for  the  refracted  longitudinal  wave  is  approached.  The  width  of  the 
peak  at  that  angle  depends  on  the  energy  at  a  b  0  and  on  the  cutoff  angle 
for  the  shear  wave.  Ttie  peak  Is  sharper  the  lower  (R/X)  at  CK  b  0  and  the 
smaller  a^.  Beyond  the  curve  dips  sharply  and  reaches  a  minimum  which 
Is  always  less  than  (R/l)^  at  a  b  0.  This  means  that  in  certain  solids  a 
shear  wave  can  be  produced  ^ich  Is  more  Intense  than  the  reflected  long¬ 
itudinal  wave  (Fig.  5)«  Whether  this  Is  possible  depends  primarily  on  the 
acoustic  Impedances  and  to  a  lesser  degree  on  the  shear  velocity. 

The  energy  of  the  shear  wave  between  a  b  0  and  Is  much 
smaller  than  Its  energy  between  and  a^,  the  maximum  depends  strongly 
on  the  values  of  and  the  width  of  the  (R/l)^  peak. 

With  this  Information,  and  since  the  two  cutoff  angles  can  be 
found  easily  from  Snell's  law  and  the  Initial  values  at  a  b  0  from  the 

reflection  coefficient,  it  Is  possible  to  estimate  the  behavior  of  the 
2  2  2 

three  curves  (R/l)  ,  (L/I)  ,  and  (S/l)  for  a  given  set  of  velocities  and 
densities  In  the  range  considered  here  without  having  uo  make  many  time- 
consuming  calculations. 


Figure  lo  Euergj  relations  et  oil-steel  bounderp  (solid  lines) 
and  water-steel  boundary  (dashed  lines). 

R  reflected  ware* 

L  refracted  longitudinal  wave. 

S  refracted  shear  wave. 


Figure  2.  Snne  as  Fig.  1  for  copper 


lO  3g>  3o 

0(-* 

Figure  3*  Sane  as  Fig.  1  for  brass. 
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Figure  1. 


Figure  2. 
Figure  3. 
Figure  4. 
Figure  5* 


Energy  relatione  at  oll-steel  boundary  (solid  lines)  and 
water-steel  boundary  (dashed  lines).  R  reflected  wave, 

L  refracted  longitudinal  wave,  S  refracted  shear  wave. 
Saote  as  Fig.  1  for  copper. 

Same  as  Fig.  1  for  brass. 

Same  as  Fig.  1  for  aluminum. 

Same  as  Fig.  1  for  magnesium. 
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